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Abstract 

We analyse the effect of the canonical quantization of the rotational mode of the 0(3) 
cr-model which includes the Skyrme term. Numerical evidence is presented that the quantum 
correction to the mass of the rotationally-invariant charge n = 1, 2 configurations may stabilize 
the solution even in the limit of vanishing potential. The corresponding range of values of the 
parameters is discussed. 



1 Introduction 

(N 

In this paper we are concerned with canonical quantization of the soliton solutions of the modified 
On ■ version of the d = 2 + 1 0(3) cr-model which includes the Skyrme term [TJ [2] (so-called baby 

Skyrme model) . The model can be viewed as a lower dimensional simplified analogue of the usual 
Skyrme model [5] , however it has its own physical application in condensed matter physics where 
the soliton solutions of the model describe spin textures in the ferromagnetic quantum Hall system 
(see, e.g., 01(7]). 

Unlike the usual Skyrme model in d = 3 + 1, a potential term is usually added to the classical 
baby Skyrme model to ensure stability of the skyrmions in two-dimensional space by the Derrick 
scaling scaling requirements. On the other hand, in the two-dimensional ferromagnetic system 
such a term describes the coupling with an external magnetic field. 

In Skyrme's original model which is considered as an effective theory of low-energy QCD in the 
limit of large number of quark colours [8j , this potential term is optional, it corresponds to the no- 
zero pion mass although its appearance might dramatically affect the structure of the solutions [9] . 
The form of the potential term in the baby Skyrme model is largely arbitrary, there are different 
families of possible models, e.g., holomorphic model which yields a simple analytic solution written 
in terms of the holomorphic functions [lOj , or double vacuum model which gives rise to the circular 
multiskyrmions [11] . Actually the choice of potential term of the baby Skyrme model dictates the 
structure of the multi-skyrmions [TH [13] . Note that apart inclusion of the potential term, there 
are other possibilities to ensure stability of the soliton solutions, for example by coupling of the 
0(3) sigma- model to a massive vector meson field [13] . 

Standard approach to the quantization of the Skyrmions involves the zero-mode quantization 
of the configuration as a rigid body [151 116j . Furthermore, the classical Skyrmion is quantized 
within the Bohr-Sommerfeld framework by requiring the angular momentum to be quantized, i.e., 
the quantum excitations correspond to a spinning Skyrmion with a particular rotation frequency. 
In the recent paper [17j an axially symmetric ansatz was used to allow the spinning Skyrmion 
to deform, it was shown that the standard Skyrme parameters are simply an artifact of the rigid 



1 



body approximation. The rotating solutions of the baby Skyrme model are also known [181 119] , 
in this case the rotationally symmetric configuration is stabilized due to radiative effects. 

Furthermore, it was suggested to treat the Skyrme model quantum mechanically, i.e., apply 
canonical quantization of the collective coordinates of the soliton solution to take into account 
quantum mass corrections [20l [21l [22l [23] . It turns out the correction decreases the mass of the 
configuration, so one can expect similar effect in the two-dimensional baby Skyrme model. 

In this paper we observe numerically that the canonical quantization of the rotational degrees 
of freedom of the baby Skyrme model produces a stable soliton solution even in the limit of 
vanishing potential term. The stabilization is achieved by including of the corresponding quantum 
corrections which appear when the canonical commutation relations are taken into account and 
decrease the mass of the soliton. 



2 The baby Skyrme model 

Let us begin with a brief review of the Skyrme model in 2+1 dimensions. In terms of the order 
parameter triplet <p — (<Pii 4>2, </>3), of the nonlinear 0(3) sigma model field subject to \4> a \ 2 = 1, 
the Lagrangian density of the baby Skyrme model has the form 

c = ^-d^ ■ d»<f> - ^(d^ X d^f - M 2 (l - 3 ). (i) 

The condition of finitcncss of the energy requires the potential term to vanish at infinity. 
This implies that on the spacial boundary <j> — > (0,0,1). Therefore the physical space R 2 is 
compactified to S 2 and such a one-point compactification ensures the topologically non-trivial 
mapping : S 2 — > S 2 . The corresponding topological charge is the homotopy invariant 

n = £ abc^ J ^X^ad^^Cpc (2) 

R 2 

There are three free parameters in the model. The constants a 2 , 1/k 2 , k 2 x [length}^ 2 and 
/i 2 x [length] 2 have the dimension of energy. On the classical level the energy scale provided by the 
parameter a 2 can be fixed in geometric units by setting a = 1, furthermore, in that case the length 
scale also can be fixed by setting k 2 to 1 but the rescaled classical parameter p remains free [T5] . 
In natural units where the Planck constant is normalized to unity we can introduce dimensionless 
radial coordinate p = —r. Then the energy functional of the static field configurations can be 
conveniently written in polar coordinates r, 9 or in dimensionless coordinates p, 9 as 



E = J rdrd9 (SL [5 r . d r <f> + ^dgcp ■ d e cj>} + ^(d r <t> x d { 



^Y + p 2 {l-h)\ (3) 



+ ^(1 -*»))■ 141 



respectively. However, when the quantum corrections to the soliton mass are taken into account, 
the energy scale is changing. Thus we cannot set a 2 = n 2 = 1 anymore and, as we will see both 
parameters are significant in that case. 

In this paper we only concerned with rotationally-invariant fields, so the corresponding 0(2) 
symmetry of the system is commonly used to re-express the field (p m terms of a profile function 
f(r) or f(p), and a polar angle 9 as 

4> = (sin f(r) cos n0, sin/(r) smn9, cos/(r)^, (5) 

where we set the global phase to be zero. f(r) is the real profile function which satisfies certain 
boundary conditions. Here we take /(0) = tt and /(oo) = 0. The integer n is actually the 
topological charge of the configuration, as one can see substituting the ansatz (O into the definition 
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(Pj). For the topological sectors with n = 1 and n = 2 this parametrization provides remarkably 
accurate approximation to the exact numerical solution of the model ([lj [2 . For higher values of 
the topological charges n > 3 this approach yields unstable circular multisoliton configurations, so 
the ground state solutions are not rotationally-symmetric although the structure of the solution 
depends on the explicit form of the potential of the model [12l [T3l [24] 

In order to apply the standard canonical quantization procedure it is convenient to re-express 
the Lagrangian {]]) in terms of the SU (2)-valued hermitian matrix fields U = <p ■ t, where r = 
(n, T2, T3) is the triplet of usual Pauli matrices 

U(fM0\ - ( cos/ ( r ) sin/We- 1 ™ 9 \ 

The asymptotic value of the field U(f(r), 9) has to tend to the matrix T3. 

Then the Lagrangian |T]) can be represented in the form similar to the usual structure of the 
Skyrme model in d = 3 + 1 

22 2 
£ = y Trd k Ud k U + ^ Tr[d k U, d t U] [d k U, d l U] - ^ Tr(l - t 3 U), (7) 

Substituting the ansatz ^ into the classical Lagrangian density ([7]) gives 

-£ ci = y /' 2 + ^|^(a 2 + « 2 /' 2 ) + M 2 (l ~ cos/). (8) 
The Lagrange density ((8]) leads to the classical Euler-Lagrange equation 

K 2 n 2 sin 2 /\ ,/ 2 , K 2 n 2 sin /cos / K 2 n 2 sin 2 / 



r(^ + ^^) + /'(« 2 + /' 



a 2 n 2 sin/cos/ 2 . 
// r sm / = U, 



(9) 



which can be integrated numerically subject of the boundary conditions imposed [2]. The equation 
© in the dimensionless coordinates takes the form 

„/ + n^sinV\ 1 /^i 1 sin /cos/ n 2 sin 2 / \ 

P 2 ■ f / 2 2 _P do) 

n sm / cos / k u 

Z-Psm/ = 0. 

p a 

22 

Thus the stability of static soliton solutions depends on a single free dimensionless parameter K % . 



3 Quantization: Momenta of inertia 

We wish to quantize the rotational degrees of freedom of baby Skyrmion by wrapping classical 
baby Skyrmion ansatz U(f(r),9) with unitary matrices A(q(t)) depending only on time t |16j 

U(q,/,0) = A(q(i))[/(/(r),0)At(q(*)). (11) 

Here, for the sake of generality, we suppose the field © of the model is embedded into the 3- 
dimensional isospace. Then the three Euler angles are associated with collective rotational degrees 
of freedom q(t) will eventually be treated as quantum-mechanical variables. The generalized 
coordinates q(i) and velocities q(i) then satisfy the commutation relations [H] 

[q\ q b ] = -i/ o6 (q). (12) 

The explicit form of the function / ah (q) will be completely determined by canonical commutation 
relations between quantum coordinates and momenta. 
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As usual, to calculate the effective Lagrangian of the rotational zero mode we have to evaluate 
the time derivative of the matrix 



U = AUA^ - Af/AtAAt, 
taking into account the commutation relations (|12[) . Explicitly, we have 

d 



1 



where 



^A(q) = C<">( g )r«A(q) = Cf } (q)A(q)r a , 



(13) 



(14) 



(15) 



and the curly brackets in (fT4|) correspond to the anticommutator. The coefficients c[ , C^°^ are 
some functions of the group parameters whose explicit form is not relevant here. 

Then, keeping only terms proportional to the square of the angular velocity in the effective 
kinetic Lagrangian density, we get 



C q (f(r),0) = iqVfc'q"' + 



iq fc Cf^)£ Q , b (/(r),^ fc) (q)q fe ' 



where the metric of the restricted configuration space of the rotational zero modes is 

9k,k'(q,f(r),e) = C^\ q )E w {f{r),e)cf\ q ). 



(16) 



(17) 



Performing explicit summation in the circular basis r + = — -^(ti+ws), to = ^t^, t_ = -1=(ti— 1T2) 
we find 



S bb ,(f(r),9) = - 



( sin 2 /e 2i " e -sin2/£^ 2 - sin 2 / \ 
3in2/^J -2 sin 2 / sin 2/^ 



V 



K 

T 



2 - sin^ / 



2 r n -2i; 



wa.2f*-j=- sin 2 /c 







c'" 8 /' 2 sin 2/ 
V2 

c ,2 f f/2 1 n 2 sin 2 / 



V cos f.r + ,,. 

Integration of the corresponding matrix 



-2/' 2 sin 2 / c "" e - f ' 2si " 2 ^ 

" 9 /' 2 sin 2/ 



/2 1 » sin 2 / 
53 



%/2 
-2in0 



V2 



(cos 2 // 



,2 n sin 2 / 
3 



Si 



/,// 



rdrd6£ w (f(r),6) 



ar 
-00 
v ai 



(18) 



gives the explicit expressions for the baby Skyrmion's momenta of "inertia" 



a = 7T 



rsin 2 / a 2 + k 2 /' 2 dr = k 2 tt / psin 2 / 1 + f' 2 ) dp = k 2 & , 



fll= 2 



• (a 2 (sin 2 / - 2) + k 2 (/' 2 cos 2 / + ^^-] 



dr. 



(19) 



(20) 



Evidently, the integral a\ in (|20[) contains a divergent term. Physically, it means the rotations 
around corresponding axes are forbidden as expected. The moment of inertia ag of a spinning baby 
Skyrmion, however, diverges only in the limit when angular velocity of the rotation q approaches 
the value of n [T8] and generally the rotations around the 3rd axis are allowed. Restricting ourselves 
to the rotations about this axis, we fix the U(l) subgroup A(q) = exp(iqr3/2) = cos |1 + i sin f T3 
of the complete isospin rotation group. 
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Figure 1: A typical energy density distribution for winding number n — 1 (a) and n = 2 (b) 
sectors. Parameters of the plotted n = 1 solution are e 2 — 28, [i s — 0, uj 2 — and correspondingly 
e 2 = 57, /Us = 0, uj 2 = for the n = 2 sector. 



Explicit calculation then gives A^(q)( J^A(g)) = ^T3, C(g) = C"(q) = i, and only the middle 
entry of the matrix (fT8|) remains. 

We are interested in a quantum-mechanical model where the corresponding quantum momen- 
tum is conjugated to the rotational collective coordinate q and it is defined as 

P = ^=^0- (21) 

The canonical commutation relation [p, q b ] — — i the allows us to define the explicit form of the 
algebra CLl 

/°°(<?) = -■ (22) 
ao 

4 Quantum Lagrangian and equation of motion 

We are now in position to evaluate the explicit form of the quantum-mechanical Lagrangian of the 
baby Skyrme model. Simple calculation yields A^A = ^r^q + g^-l and the quantum Lagrangian 
is given by 

L q = J C q [f(r)] = £--AM (23) 

where 

» „ , 7T f , / o . 9 „ k 2 / „, 9 . 9 „ n 2 sin 2 f 2n 2 sin 4 f\\ . 

= "8^ J ^ [ a Sm / + 32 ( 32/ Sm ' + ^ + ^) J (24) 

7T /" / , 1 / „ , /9 ?i 2 sin 2 / 2n 2 sin 4 f\\ AM , x 

is the quantum mass correction which appear when the commutation relations (|12[) are taken into 
account. 

We define the angular momentum operator as 

J = -]^{PiC^ l {q)} = -a q. 



Its eigenstates are the vectors 

\u) = exp(-iwg)|0), 

where to is an integer which enumerates the irreducible representations of the U(l) group. 

The total effective Hamiltonian corresponds to the complete Lagrangian L = L c i + L q which 
includes both classical and quantum mechanical parts: 

H = \{p, q}-L = ^--L cl + AM > ( 26 ) 

In the framework of the Bohr-Sommerfeld quantization of the rotational zero mode of the baby 
Skyrmion, the quantum mass term is absent [18j . This picture corresponds to the rigid rotation of 
the Skyrmion with fixed profile of the function f(r). However, since the quantum mass correction 
AM turns out to be negative it might stabilize the baby Skyrmion solutions as it happens in the 
three-dimensional Skyrme model |25j . This is the most interesting problem to study, so henceforth 
we shall mainly consider the limiting case of vanishing potential term p — 0. 

Indeed, let us consider the equations which corresponds to the minimization of the total energy 
functional. Varying it we obtain rather cumbersome integro-diffcrcntial equation in dimcnsionless 
coordinates on the profile function f(p) replacing its classical counterpart ([9]): 

f"{p) ( 2rt2si ; 2/(P) +2p(l + Z M g S inV(p))) 
+ /' 2 (P)(^^ + P^sin2/( P ) 

+ /'(p)(-^^ + 2Z(l + ^sin 2 /(p))) ^ 
- sin 2/00 P^A + + 4 sin 2 f(p)(l - Z)) 

V p 

-2Z/x 2 psin/(p)=0 

Here the usual boundary conditions on the function f(p) arc imposed and we introduce the 
shorthand notations for the dimcnsionless quantities 

^ 2 2 

Z = !~ o8„.2..2=2 ; M' = T^F ( 28 ) 



and 



where 



(pf' 2 (p) 



2Za 2 K 2 o 2 V4 2 6 a 

n 2 sin 2 /(p) 2n 2 sin 4 /(p) 



P P 
As p — > oo, the equation (|27p reduces to the asymptotic form 



dp. (30) 



Pf"(p) + f'(p) - ™ 2 pf(p) = 0. (31) 

where the quantity m 2 = p? s + p 2 , corresponds the asymptotic mass of the excitations. Thus the 
leading term in an asymptotic expansion of the function f(p) is given by 

e -mp , i N 
M = Cl — (l-—). (32) 

and the solution of the quantum-mechanical model remains exponentially localised. The constant 
of integration C\ which appears here will be determined later from the results of the numerical 
calculations. 
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5 Numerical Results 



The integro-differential equation (j2"T|) can be solved numerically by shooting method. The initial- 
ization of the algorithm requires trial values for all the integrals ao, p 2 d and Is that appear in the 
equation to be specified. Estimates of these can be obtained by employment of the classical profile 
function f(p) of the corresponding winding number n baby Skyrmion. If we drop out the explicit 
mass term fi s = (the most intriguing case), the stable soliton solution can be obtained only 
for some range of the values of the parameter 1/e 2 = a 2 n 2 . This can easily be seen by noticing 
that in the limit e 2 — > the quantum equation (|2"T)) is formally reduced to the classical equa- 
tion dTUJ) . However in this limit the classical chiral function f(p) decays as ~ - and, consequently 
the moment of inertia given by the integral ao (fH?)) diverges [TH HH] • Physically it means that as 
e 2 decreases, the quantum baby Skyrmion slows down and the quantum correction to the soliton 
solution becomes negligible, so they cannot stabilize the configuration. 

On the other side, the limit of the large values of the parameter e 2 , i.e., the case of relatively 
large quantum corrections, implies small values of the quartic stabilizing term, so in this limit 
the soliton solution becomes unstable. Thus, we may expect existence of a window of values of 
parameter e 2 for which we may get a quantum-mechanical rotated solution with a non-vanishing 
potential term ji s ^ 0. 

The solution of integro-differential equation then proceeds as follows. Certainly, for some 
intermediate range of values of the parameter e 2 , the values of the integral ao can be obtained by 
introducing an integration cutoff. These approximated values then can be used as an input for 
the next step of numerical iteration over the entire range of values of the radial variable. 

Shooting from the point /9 max (where /(/?) assumed to be of the form (J22)), to the point p mm 
(here f(p) = f(p m i n ) — (p m in ~ p)/'(Pmin)) and varying the only unknown constant C\ in (1321) 
yields a continuous family of solutions, which satisfies the required topological boundary conditions 
/(0) = 7r and /(oo) = 0. Typically, pmax G [5, 12] and we set p m i n ~ 10~ 3 to obtain solutions with 
topological numbers n = 1, 2, 3. 

For larger values of the topological charges n > 4 the multisoliton solution profile function 
/ slowly varies in the vicinity of the origin, so the energy density distribution is getting more 
extended 0. Technically it means the value of p m i n should be increased up to lCP 1 to keep the 
numeric algorithm stable. 

Once the profile function is found, it can be used as an input for the next step of the iteration 
procedure. Then we recalculate all required integrals again and repeat the same procedure until 
all the integrals do converge to some stable values. 

The analytical and numerical calculations are performed with Mathematica [26] . Typically 
20-60 shoots are enough to achieve high precision numerical solution of the differential equation 
and 1000-5000 iterations is needed to ensure all the integrals are definitely converge to some set 
of fixed values in the case when we set p s = 0. 

In figure we presented the results of numerical calculations for some particular values of the 
parameters of the model. The profile functions and the energy density distributions of the soliton 
solutions with n = 1 ... 4 are shown in Fig. [2j Here we take the value e 2 = 30 in the sector with 
topological charge n — 1 and e 2 = 40 for the n = 2 soliton and e 2 = 80 for the solution solutions 
with n = 3, 4. 

Note that the maxima of the energy density distribution of the n = 2 . . . 4 solitons are shifted 
away from the origin. The size of the solutions increases as the winding number n increases. 

In Figures [3l|4] we present the properties of the p s = soliton solutions within the sectors n = 1 
and n — 2, respectively. Evidently, the soliton's energy density distribution becomes more and 
more spread out, i.e., the characteristic size of the soliton increases as the parameter e 2 decreases. 

1 Remind that the rotational symmetry of the classical baby Skyrmion solution holds in the sectors with n = 1, 2 
providing an absolute minimum of energy functional there. However the effect of rotation may affect the structure 
of the discrete symmetry solutions with higher values of the topological charge n > 3 increasing degree of symmetry 
up to rotationally invariant ansatz |J5J. Although we do not address this issue here, we hope to investigate this 
transition in our future work. 
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Figure 2: Typical profile functions (dashed lines ) and the distribution of the energy density (solid 
lines ) of the soliton solutions with winding numbers n = 1, . . .4 and ji s = 0, u> 2 — 0. Numerical 
values of the parameters e 2 , /i 2 ,, the integrals Ib and ao and the total integrated mass of the 
solitons are presented in the graphics legend. 




Figure 3: The profile functions (dashed lines ) and the distribution of the energy density (solid 
lines ) of the n = 1 solutions with = 0, lu 2 = 0. Numerical values of the parameters e 2 , /it 2 ,, 
the integrals Ib and ao and the total integrated mass of the solitons are presented in the graphics 
legend. 
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Figure 4: The profile functions (dashed lines ) and the distribution of the energy density (solid 
lines ) of the n — 2 solutions with fi s = 0, uj 2 = 0. Numerical values of the parameters e 2 , /it 2 ,, 
the integrals Ib and ao and the total integrated mass of the solitons are presented in the graphics 
legend. 



The lines marked with 1 in both figures represent the solution which corresponds to the maximal 
possible value of the dimensionless parameter e 2 . Our analysis shows that a stable solution does 
not exists when this parameter increases beyond this critical value, i.e., quantum correction to 
the soliton mass becomes to small to stabilize the solution in the limit of the vanishing potential 
term. Evidently, this critical value increases if /i s is taking to be non-zero. 

On the other hand, increasing of the parameter e 2 results in increasing of the values of the 
momenta of inertia of the quantum soliton. In other words, the quantum correction to the soliton 
mass become large and the characteristic size of the configuration decreases. However, as we 
can see from plots presented in Figs. [3j IH further increasing of the values of the parameter e 2 
yields negative values of the energy density distribution (see the curves marked with 4), so the 
corresponding solutions should be considered rather an artifact of the numerical calculations. 

Therefore if fi s = 0, the physical soliton exists for some domain of values of the parameter e 2 
which is restricted from both sides. Evidently, the quantum correction to the mass of the soliton 
cannot be very large, it has to be of order of about a few percent of the classical mass. The upper 
bound on the range of values of the parameter e 2 corresponds to the ratio AM/M c i ass ~ 0.4 which 
seems to be too large value for a quantum correction. 

Conclusion 

We investigate the effect of quantization of the rotating baby Skyrmions beyond the usual Bohr- 
Sommerfeld framework and the rigid body approximation. Our results indicate that the canonical 
quantization of the rotational collective coordinate of the model may stabilize the soliton solution 
even in the limit of vanishing potential term. We have checked that the corresponding radially- 
symmetric solutions exist in the sectors with winding numbers n — 1 , ... 4 for a range of values 
of the parameters of the model. The energy density distribution of the quantum baby Skyrmion 
is exponentially localised as we have seen from the asymptotic formula (|32| . The shape of the 
quantum soliton depends on the value of the quantum correction to its mass. At the critical values 
of the parameters of the model the configuration becomes unstable or the quantum correction to 
the mass is too large and the corresponding energy distribution becomes negative. 
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It would be interesting to investigate the structure of the quantum solitons in different topo- 
logical sectors with various potentials included. In particular, a question arises about the shape 
of the minimal energy solution. On the one hand, it is known the choice of the potential strongly 
affects the structure of the multisoliton solutions of the model [TU Q2J EH HI] • On the other hand, 
rotation of the field configuration may result in the restoration of the radial symmetry, so this 
problem is worth investigating further. 

This work is supported by the Science Foundation of Ireland in the framework of the Science 
Foundation Ireland (SFI) Research Frontiers Programme (RFP) project RFP07/FPHY330. We 
are grateful Tigran Tchrakian for valuable discussions and remarks. 
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